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Abstract 

We analyze the constraints placed on individual, flavor diagonal CP-violating phases in the minimal 
supersymmetric extension of the Standard Model (MSSM) by current experimental bounds on the 
electric dipole moments (EDMs) of the neutron. Thallium, and Mercury atoms. We identify the 
four CP-violating phases that are individually highly constrained by current EDM bounds, and 
we explore how these phases and correlations among them are constrained by current EDM limits. 
We also analyze the prospective implications of the next generation of EDM experiments. We 
point out that all other CP-violating phases in the MSSM are not nearly as tightly constrained 
by limits on the size of EDMs. We emphasize that a rich set of phenomenological consequences is 
potentially associated with these generically large EDM-allowed phases, ranging from B physics, 
electroweak baryogenesis, and signals of CP-violation at the CERN Large Hadron Collider and at 
future linear colliders. Our numerical study takes into account the complete set of contributions 
from one- and two-loop EDMs of the electron and quarks, one- and two-loop Chromo-EDMs of 
quarks, the Weinberg 3-gluon operator, and dominant 4-fermion CP-odd operator contributions, 
including contributions which are both included and not included yet in the CPsuperH2.0 package. 
We also introduce an open-source numerical package, 2LEDM, which provides the complete set of 
two-loop electroweak diagrams contributing to the electric dipole moments of leptons and quarks. 
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I. INTRODUCTION 



The minimal supersymmetric extension of the standard model (MSSM) features a large 
number of additional CP-violating sources compared to the SM, yielding an extremely rich 
array of possible phenomenological consequences [H, H]- To mention a few, the new sources 
of CP violation can participate in the generation of the baryon asymmetry of the universe in 
the context of supersymmetric electroweak baryogenesis They can generate observable 

CP asymmetries at the CERN Large Hadron Collider (LHC) jsj-S] as well as at future linear 
colliders js-l^l. They could manifest themselves in B physics (see Ref. [l3-[l6| and the recent 



results from the Tevatron DO Collaboration [iTj] and CDF Collaboration[l8|). In general, 



CP-violating phases induce rather large contributions, compared to current experimental 
sensitivity, to the electric dipole moment (EDM) of the neutron and of atoms [L^. A 
large CP-violating phase could also have an effect on many CP-conserving observables such 
as the mass spectrum, the production rate and decay branching ratios of SUSY particles, 
especially the lightest neutralino as a dark matter particle (including direct and indirect 
detection rates, and the relic density) (2l|-23| and Supersymmetric Higgs scalars^, 25 . 



All these potentially interesting observable signatures are however highly constrained 
by the fact that no permanent EDM has ever been experimentally observed. The current 
most stringent bounds are on the EDMs of the neutron [2^, and on that of Thallium 
[2?! and Mercury [2^ atoms^. Specifically, the current experimental constraints on the 
aforementioned quantities read: 

|rfn| < 2.9 X IQ-^^e cm (90%C.L.), 
I^Til < 9.0 X IQ-^Sg cm (90%C.L.), 

I^Hgl < 3.1 X IQ-^^e cm (95%C.L.). (1) 

Signatures of CP-violation compatible with the bounds described above have so far been 
discussed in the context of scenarios where one avoids EDM bounds and keeps relatively large 
phases with light super-particle masses via one or more cancelations among the various terms 
contributing to the EDMs (I, 29-35|. This "cancellation" scenario occurs typically in some 



very fine-tuned region^ of the MSSM parameter space. The goal of the present analysis is 
instead to examine in detail, beyond the possibility of cancelations between different CP- 
violating sources^, how each individual phase is constrained by current EDM bounds and to 
study under which conditions large phases are phenomenologically allowed while also keeping 
the relevant mass scales relatively light and thus phenomenologically interesting. Instead of 
imposing any universality condition as in the supergravity (SUGRA) model, we entertain 
here the possibility of general non-universal soft terms (as is the case in many string-inspired 
models, see e.g. |37|), and thus keep all the soft terms in MSSM as independent variables. 
For concreteness, we focus on flavor diagonal phases. We defer a study of EDM bounds on 
flavor non-diagonal phases to future work. One of our main results is that among all MSSM 
CP-violating phases, only four are in fact strongly constrained by current EDM bounds, 



^ For recent reviews on EIQM searches and their imphcations for the MSSM, see, e.g. Refs. 0, 0, 0|- 



^ It is estimated in Ref. 



that a minimal degree of fine-tuning at the level of 10 ^ is needed. 



There are still cancelations among different contributions to EDMs from the same source in some region 
of parameter space. We emphasize that our main focus is on the general trend of constraints without 
cancelation, and always point it out when we encounter such cancelation region. 
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and we study the correlations among them. We comment on the possible phenomenological 
implications of those phases that are not strongly constrained. 

Part of our numerical study is based on the CPsuperH2.0 package which includes 

the complete contributions from one-loop supersymmetric contributions to the EDMs of 



33, 40 



the electron and of quarks, one-loop and two-loop Chromo-EDMs of quarks [31 
6-dimensional 3-gluon Weinberg operator j4l[, and dominant contributions to 4-fermion 
CP-odd operators [4^, but only a subset of the contributions of two- loop EDMs of electron 
and quarks j4o|. There exist additional Higgs exchange-mediated chargino-neutralino 2- 
loop contributions to electron and quark EDMs that can become dominant in the limit 
of heavy sfermions [43!] and that are not included in CPsuperH2.0 package. We here take 
into account these contributions as well, therefore featuring the dominant 1-loop and 2-loop 
contributions to all the 6-dimensional CP-odd operators that generate EDMs of the neutron 
and of Thallium and Mercury atoms^. The numerical code where all these new contributions 
are collected is called 2LEDM, and is currently available from the authors upon request. The 
2LEDM code currently includes an interface to the FeynHiggs^ package (version 2.6.5) ji^ ]. 
We plan in the near future to set up a webpage for easier download of 2LEDM and of related 
tutorials and to include an interface to the CPsuperH2.0 package. 

There remains an order one theoretical uncertainty with neutron and Mercury EDMs. 
For the neutron EDM, the uncertainty arises from hadronic physics, while for the Mercury 
EDM the source of uncertainty is associated with (i) atomic physics in extracting the nuclear 
Schiff moment from c^Hg, (ii) the nuclear physics going into extracting T- and P-odd pion- 
nucleon couplings q-^nn from the Schiff moment, and (iii) the hadronic physics in computing 
the Qt^nn in terms of quark Chromo-EDM operator, Weinberg three-gluon operator, and 
CP-violating four fermion operators^. In utilizing the CPsuperH2.0 package to estimate 
the relevant EDMs from the EDMs of quarks and leptons we are relying on QCD sum rule 
computations 42|, |45|, |46| of strong interaction matrix elements. For a discussion of the 
systematic uncertainties in our results for these quantities when different hadronic model 
approximations are employed we refere the reader to e.g. 



Ref. 29 



We note also that 

the CPsuperH2.0 code relies on the computations of the nuclear Schiff moment reported 



in Ref. |47[ and does not take into account the recent computations of Refs. |48 



r 

Mi- The 

i99Hg to 



latter two computations give an enhanced sensitivity of the nuclear Schiff moment of 
the isoscalar, T- and P-odd pion-nucleon coupling, g^^ as compared to Ref. 47|. Moreover, 



Ref. 49 1 finds that the sensitivity to (7!^^ may be reduced by a factor of three to ten, depending 
on the type of interaction used. 

In general, g^^ are dominated by quark Chromo-EDMs rather than the Weinberg three 
gluon operator (the contribution is suppressed by mg), while in the MSSM with large tan/3, 
contributions from the four fermion operators may be important 19|, |20|. The QCD sum rule 
analysis implies that (7!°^ is five times less sensitive to the sum of up and down quark Chromo- 
EDMs than g'^^ is to their difference. Thus, we would only expect the stronger sensitivity of 
the Mercury Schiff moment to g'^^ to be important in small corners of the MSSM parameter 
space where the difference of quark Chromo-EDMs is highly suppressed compared to their 
average value. The possible suppression in sensitivity to g'^^ is a potentially more serious 



^ The four fermion operators in the MSSM are technically dimension eight. 
^ http : //www ■ f eynhiggs . de/|| 

^ In the MSSM, the Chromo-EDM operator typically gives the dominant contribution to the g-^NN ■ 
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issue. Thus, one may need to relax the constraints we obtain on CP- violating phases that are 
driven by the ^^^Hg results in light of on-going theoretical nuclear structure developments. 
With these caveats in mind, we summarize our main findings here: 

(a) A primary impact of the new ^^^Hg result is to impose significantly more stringent 
constraints on the relative phase 03 between the gluino soft supersymmetric-breaking 
mass and the /i parameter (see below), while generating a strong correlation between 
this phase and the phase of the soft-breaking triscalar couplings involving first gener- 
ation sfermions. 

(b) The neutron and Thallium EDM limits have a stronger impact on the relative phase 
02 between the wino soft mass parameter and fi than does the ^^^Hg bound, but at 
present there does not exist any strong correlation between 02 and other phases. 

(c) A future neutron EDM limit that is roughly 100 times stronger than present would both 
tighten the present correlations between 03 and the triscalar phases while inducing 
strong correlations between 02 and other phases. 

(d) In the limit of heavy first and second generation sfermions, the "bino" phase 0i is 
essentially unconstrained by present EDM bounds. A future neutron or electron EDM 
measurement with ~ 100 times better sensitivity would probe the impact of this phase 
at a level of interest for cosmology. 

In the remainder of the paper, we organize the discussion of our analysis leading to the 
findings above as follows: In section [TTl we give a general discussion about the CP- violating 
phase structure of the MSSM, we address how each phase impacts the various EDMs, and 
we outline the eneral setup of our analysis. In section IIIH we investigate in detail how each 
phase is constrained by current EDM bounds, we study the correlations between the various 
EDM bounds on the most strongly constrained phases, namely 02, 03, and (f)u,d, and we 
discuss the phenomenology implications of the other loosely constrained phases. Finally, we 
devote section IIVI to our summary and conclusions. 

II. CP- VIOLATING PHASES IN MSSM AND THE SETUP FOR ANALYSIS 

The Minimal Supersymmetric Extension to the Standard Model of particle physics intro- 
duces a plethora of new and unknown parameters. Many of these parameters are connected 

TABLE I: Summary of how the CP-violating sources in MSSM generate various CP-odd operators 
at one-loop and two-loop level. 



CP-violating phases 


one-loop contribution 


two-loop contribution 


0e,M,d 


il—loop 7l~loop 


no 


0/i,c,s 


no 


no 


4>T,t,b 


no 




01,2 






03 


jl—loop 71— loop ^ 

%4 ' 
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to new sources of CP or flavor violation, or both. Although EDMs could, in principle, be 
induced by all CP-violating parameters including both flavor-conserving and flavor- violating 
ones, they are most sensitive to flavor-conserving CP- violating phases, including those asso- 
ciated with the bilinear coupling b and Higgsino mass term /x in the Higgs-Higgsino sector, 
the soft-supersymmetry breaking Majorana masses Mi, M2, and M3 in the gaugino sector, 
and the trilinear couplings ^/ in the sfermion sector. 

This notwithstanding, not all the new CP violating phases appearing in the MSSM are 
physical. In fact, there exist two transformations that can be employed to rotate away two 
phases [soj. We choose a convention where fi and b are real, and the remaining phases 
mentioned above are all physical. In particular, the physical phases include the phases 
4^1,2,3 of the gaugino masses Mi^2,3, and the phases (pu,d,e, (pc,s,^i, and 4>t,b,T of the sfermion 
trilinear couplings Au^d,e, Ac,s,^J., and At^b,T, respectively. As shown in Table HI these phases 
play different roles in generating various CP-odd operators, including the electron EDM 
de, quark EDM dg and Chromo-EDM dq, the Weinberg 3-gluon operator d^'^, and the 4- 
fermion CP-odd operator Cff/'^. These CP-odd operators are responsible for the EDMs of 
the neutron, as well as of that of the Thallium and Mercury atoms, as summarized in Table 
im In particular, in the MSSM the Thallium EDM is dominated by the electron EDM 



operator de, and possibly by the four-fermion operator Cf f/ if tan/3 > 30 [42[; the neutron 
EDM, which we compute here using QCD sum rule results |45| , mainly stems from the EDM 
and chromo-EDM operators of the u and d quarks, du^d and du^, and from the 3-gluon term 
(i^*^; lastly, the Mercury EDM is generated primarily by the chromo-EDM operators du^d 
[l9| . A combination of Table H] and Table [III provides information on how each CP- violating 
phase is constrained by which experimental EDM bound. 

Among all contributions, some of the dominant ones stem from the one-loop induced 
EDM and Chromo-EDM operators de, d^^di and d^^d- These contributions always involve the 
first-two generations of sleptons and squarks, and therefore are asymptotically suppressed 
in the limit where these scalar fermions are very heavy [5l| . Obviously, the effect of the CP- 
violating phases (pu, (pd, 4>s, (pc, (pe, <Pfj, from the first-two generations sfermions are completely 
suppressed in this situation, and would not show up in any other observable signature. In 
contrast, the effects of other phases, including 0i,2,3 in the gaugino sector and (pt, (pb, (pr 
in the third-generation sfermion sector (thanks to larger Yukawa couplings), are not as 
strongly suppressed in the decoupling limit of heavy first and second generation sfermions, 
and they might induce interesting effects that could manifest themselves at colliders or in 
other experiments sensitive to CP violation. 

With these considerations in mind, we study cases where the sfermion masses for the first 
two generations are either light or heavy. For each case, we explore in detail the mass- and 
tan^-dependence of the EDM bounds on each individual phase. We choose a set of mass 

TABLE II: Summary of relevant CP-odd operators of neutron, Thallium, and Mercury atom EDMs. 



dn 


dTi 


dug 


du,d, du,d, d^'^, Cff 


de, Cff 


de, du,d, Cff 



For the specific form of each of these operators, see, e.g. Refs. [2 
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parameters corresponding to a light spectrum as the reference point (we call this "Case 
I"); we then study the effect on the constraints as the relevant mass scales increase (the 
limiting case is indicated as "Case II"). The reference values we choose for the relevant 
supersymmetric parameters Mi 2,3, fi, Aj oi all flavors /, the charged Higgs mass Mh±, and 
the third generation sfermion masses mi^^R-i in both cases are as follows: 

|Mi| = 150 GeV, IM2I = 250 GeV, \M^\ = 550 GeV, 
= 225 GeV, \Af \ = 175 GeV, M^i = 500GeV, 
= mns = 200 GeV, (2) 

We set the first-two generation sfermion masses, in the two cases, to: 

CASE I : 2 = niR^ ^ = 200 GeV, (3) 
CASE II : m^, 2 = ^Ri,2 = 10 TeV. (4) 

In the study of each phase, we look at how EDM constraints are affected by changes in the 
relevant mass scale, keeping all other masses set to their reference values. 

III. DETAILED ANALYSIS 

In the present study, we classify the MSSM CP-violating phases under consideration into 
two groups, 

(i) phases in the higgsino-gaugino sector: 0i, 02, and ^3, and 

(ii) phases in the sfermion sector: (pa, 4>e, (pc, (ps, 4>fj.^ 4>t, (pb, 4>r, 
and we study the two groups of phases individually. 



A. Phases in the higgsino-gaugino sector: (pi, (p2i and cj)^ 



These three phases contribute to the EDM and Chromo-EDM operators both at the one- 
and two-loop level; we therefore discuss here the constraints on them for both cases: (I) 
with light and (II) with heavy first-two sfermion generations. 



Light sfermions 



In the case where the first-two generation sfermions are light and one-loop EDMs and 
Chromo-EDMs are not suppressed, the phases (pi and 02 induce contributions to df., while 
01,2,3 generate {du,d, du^d), and the four-fermion CP-odd operators Cff. In addition, 03 also 



induces a non-zero contribution to the 3-gluon operator d^'^ at two-loop order 53 



We find that as far as the constraints on 0i are concerned, the Mercury EDM, with its 
newest experimental bound, puts much more stringent bounds - namely by a factor of 10 
or more on most of the parameter space we consider - than those from the neutron and 
Thallium EDM limits. We illustrate the resulting limits in Figure [H top panels. For the 
phase 02, instead, the Mercury EDM puts less stringent limits, by a factor of a few, than 
the current Thallium and neutron EDM bounds, illustrated in Figure [H middle panels. We 
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FIG. 1: Constraints on the CP violating phases (pi, (t)2 and </»3 (from top to bottom) versus {mL,R)i 
(left panels) and versus Mi^ 2, 3 (right panels) from experimental limits on the Mercury (black 
solid lines), Tallium (red dotted) and neutron (blue dashed) EDMs, in the case of light first- two 
generations of sfermions (case /, see Eq. [3]). 
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FIG. 2: Curves of constant values for the Thallium (red) and neutron (blue) EDM as a function of 
Ml and Mjj± in the case of heavy first- two generations of sfermions {case IT) as in Eq. HI Because 
current EDM limits do not constrain (pi in this case, curves correspond to representative future 
EDM sensitivities. Solid and dashed curves correspond, respectively, to tan/3 = 3 and 60. 
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FIG. 3: The M2- and Mjj±- dependent constraints on (/>2 from neutron and Thallium EDMs in the 
case of heavy first- two generations of sfermions (case IT) as in Eq. HI 



find that this is due to the cancellation of 02 contributions to the electron EDM and quark 
chromo-EDMs that generate Mercury EDM. In fact, in a different mass region where the 
cancelation is not significant, the Mercury EDM constraint on 02 may be comparable or even 
stronger than the current Thallium and neutron EDM bounds. For 03, the new Mercury 
EDM bound puts a stronger constraint than the neutron EDM bound, while the current 
Thallium EDM bound is not stringent enough to put any constraint on 03, due to suppressed 
contributions from Cff (see Figure [H, lower panels). 

Our results as a function of the relevant mass scales are summarized in Figure [H In the 
upper penals, we show how the Mercury EDM constraint on 0i depends on tan/3 (set to 
3 and to 60 in the upper and lower curves, respectively) and on the relevant mass scales 
{rriL^Rji (left panel) and Mi (right panel). The most important impact of 0i on the Mercury 
EDM is through the neutralino one-loop contribution to the quark Chromo-EDM ^, where 
the external gluon is only attached to squarks in the loop. In this case, the dependence on 
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FIG. 4: The {mL^ji)^- and M3- dependent constraints on (/>3 from neutron EDMs in the case of 
heavy first- two generations of sfermions (case IT) as in Eq. HI 



the sfermion and gaugino masses are somewhat non-trivial: The allowed values for 0i grow 
quickly with increasing {mL,R)i (left panel), but rather slowly with increasing Mi (right 
panel). As {mL^R)i increases to 1.2 TeV, 0i can be as large as 7r/2(0.037r) for tan/3 = 
3(60), indicating that a larger value for 0i is allowed for smaller tan/3 when the one- loop 
contributions dominate. 

In the two middle panels of Fig. [H we show how the neutron (blue dashed lines) and 
the Thallium (red dotted lines) EDM constraints on 02 depend on tan/3 and on the relevant 
mass scales {1711^^)1 and M2. We notice that the current neutron and Thallium EDM bounds 
put comparable constraints on the wino phase 02 • The most important contributions from 

02 to the neutron and Thallium EDMs are through the chargino one-loop contribution to 

± 

quark and electron EDMs d^^^, where the external photon is attached to both sfermions 
and charginos in the loop. In this case, the dependence on the sfermion and gaugino masses 
resemble the case of the neutralino loop. A comparison of the left and right panels shows that 
the allowed values of 02 grows a little faster with increasing {mLj^)i than with increasing 
M2. For tan/3 = 3(60), the phase 02 is allowed to be within 10~^7r(0.6 x 10~^7r) for (m^^ /j)i 
up to 1.2 TeV, and 0.5 x 10-^(0.3 x 10"%) for M2 up to 1.25 TeV. 

In the lower panels of Fig. [H we finally show how the constraints on 03 depend on tan/3 
and relevant mass scales {mL,R)i (left) and M3 (right). The black line indicates the Mercury 
EDM constraint, while the blue dashed line refers to the neutron EDM. We observe that 
the constraints on 03 do not monotonically increase with mass scales, due to non-trivial 
cancelations among du,d, du,d, and d^^ in their contributions to the neutron EDM, and to 
cancellations among different contributions to d^^d, which dominate the Mercury EDM. How- 
ever, there is no common region where both the neutron and the Mercury EDM constraints 
are suppressed; hence, the cancellations among various 03-dependent contributions never 
open a region in this portion of parameter space region where this CP-violating phase can 
be large, independent of the values of the other phases. For tan/3 = 3(60), the phase 03 is 
allowed to be within 10"V(0.5 x IO^^tt) for (771^,^)1 up to 1.2 TeV, and 10~V(0.2 x lO^^vr) 
for M3 up to 1.55 TeV. 
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Heavy sfermions 



In the case where the first-two generations of sfermions are heavy (case II), and therefore 
the one-loop EDMs and Chromo-EDMs are suppressed, the gaugino phases 0i and 02 still 
induce a non- vanishing 4-fermion CP-odd operator Cff due to loop-induced mixing between 
the CP-even and CP-odd Higgses that are exchanged between the fermions. In addition, 0i 
and 02 induce de and du^d at the two-loop level, while 03 induces 3-gluon operator d^'^ and 
Chromo-EDM du,d at the two-loop level. 

Fig. [2] shows that in case II (heavy sfermions) the phase 0i is not constrained by cur- 
rent EDM bounds at all. This has important consequences for scenarios where the baryon 
asymmetry in the universe is generated via the mechanism of electroweak baryogenesis: a 
non-vanishing and large enough 0i can generate the observed baryon asymmetry via the 
bino-driven scenario in the context of electroweak baryogenesis, even if 02 = 8. In Fig. H 
where we set 02 = 0, we show contours of constant values for the neutron and Thallium 
EDMs that are below the current bounds, at tan /3 = 3 and 60. Notice that, in contrast 
to the situation for tan /3 = 60, at tan /3 = 3 cancellations occur between the WW and 



the W contributions, the only non- vanishing graphs for 02 = [43[. In particular, 
Ref. pointed out that these two contributions have opposite signs, and if tan/3 has a 
value such that the two are comparable, then cancellations are possible, and indeed they 
occur for Mi ~ 320 GeV and for ttih ~ 420 GeV in our setup, as shown in Fig. [2l Similar 
cancellations do not take place at large tan P, where the W^H^ contribution dominates. 
If instead the wino phase 02 is non- vanishing (Fig. 12]), two-loop contributions are much 



larger (see Ref. j43| for a discussion on what makes wino-driven two-loop EDMs typically a 



factor of 50-100 larger than bino-driven ones). Specifically, with our choice of parameters, the 
limits on 02 as a function of M2 from the EDM bounds of Thallium atom and neutron range 
from 0.006 tt for small M2 = 200 GeV up to 0.02 tt for M2 ~ 1 TeV. A similar dependence is 
found for the second mass scale entering the two-loop contribution, namely the heavy Higgs 
sector, where we find that 02 must be smaller than 0.007 vr for small mH± ~ 500 GeV and 
than 0.03 vr for m^i ~ 1500 GeV. 

In case II, the phase 03 induces EDMs of neutron. Thallium, and Mercury. Since it does 
not induce an electron EDM d^, its contribution to Thallium EDMJs highly suppressed. The 
phase 03 could induce sizable Mercury EDM through generating du^d at two-loop. However, 
the most stringent constraint comes from the current neutron EDM bound. In Fig. HJ we 
show how the constraint on 03 from neutron EDM bound depends on tan/3 and relevant 
mass scales {mi^Rjs and M3. We observe that for tan/3 = 3(60), the 03 is allowed to be as 
large as 7r/2(0.027r) for {mL,R)i as heavy as 1.2 TeV, and 0.3 x 7r(0.027r) for M3 < 1.55 TeV. 



B. Phases in the Sfermion sector: 0e, 0^, 0^, 0^, 0C5 0S5 0t) 06) 0r 

We discuss here how CP violating phases in the sfermion sector are constrained generation 
by generation. The phases associated with the first sfermion generation, namely 0e, 4>u and 
0d, yield significant contributions to the electron and u- and d-quark EDMs, as well as to 
the Chromo-EDM, all only at the one-loop level. Thus, we only discuss here the case where 
the first-two generations of sfermions are light (case I) as in the case where sfermions are 
heavy all of these contributions asymptotically vanish. 

Besides contributions to the 4-fermion CP-odd operators Cff which are in general small, 
the phase 0e contributes to the electron EDM, and is therefore constrained by the Thallium 
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FIG. 5: The {mi^^i- and Mi- dependent constraints on (j)^ from Thallium EDM in the case with 
light first-two generations of sfermions as in Eq. [3J The tan/3 = 3 is used in making the plot. The 
constraints corresponding to other values of tan/3 are not shown, as the tan/3 dependence is found 
to be rather weak. 




FIG. 6: The {rriL^jiji- and M3- dependent constraints on (pu from neutron and Mercury EDMs 
in the case with light first-two generations of sfermions as in Eq. [3l The tan/3 = 3 is used in 
making the plot. The constraints corresponding to other values of tan/3 are not shown, as the tan/3 
dependence is found to be rather weak. 



EDM. The phases 0„ and 0^ contribute to the EDM and to the Chromo-EDM of quarks, and 
are therefore constrained by both the neutron and the Mercury EDMs. The relevant mass 
scales in these one- loop contributions are (^l,/?)! and Mi 2,3- However, since the one-loop 
contribution involving neutralinos dominates d^, the constraint on 0e depends mainly on 
{^l,r)i and Mi. On the other hand, constraints on the squark CP violating phases (pu and 
(pd mainly depend on {mL ji)i and M3, since here the one- loop contribution involving gluinos 
dominates du^d and du^- The dependence on tan/3 is found to be rather weak. 

As shown in Fig. O the allowed value of 0e grows much faster with increasing {nii fiji 



11 




{mL,R)i [GeV] 



M3 [GeV] 



FIG. 7: The {mi^fiji- and M3- dependent constraints on (pd from neutron and Mercury EDMs 
in the case with hght first-two generations of sfermions as in Eq. [31 The tan/3 = 3 is used in 
making the plot. The constraints corresponding to other values of tan/3 are not shown, as the tan/3 
dependence is found to be rather weak. 
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FIG. 8: The (mL,_R)3-dependent constraints on (j)t from neutron, Thallium, and Mercury EDMs. 



than Ml, due to the dominant neutralino one- loop contribution. Constraints on this phase 
are rather weak and they reach 7r/2 for wil,_Ri of only 500 GeV, and for Mi of 1.15 TeV. 

As far as constraints on (pu and (f)d are concerned, the new experimental limit on the 
Mercury EDM places a more stringent constraint (by a factor of a few) than the current 
neutron EDM bound. As shown in Figs. [6] and [71 which display the {mL,R)i- and Ma- 
dependent constraints on and 0^, respectively, the phase (pui't'd) is constrained to be 
< 0.077r(0.047r) for {mL,R)i < 1.2 TeV, and 0.l7r(0.027r) for M3 < 1.55 TeV. 

The phases associated with the third generation sfermions 06, and 0,- do not directly 
induce an electron or quark EDM, or a Chromo-EDM at the one-loop level, and therefore 
there is no difference here between the case with light and heavy first-two generations of 
sfermions. The third generation CP violating squark phases contribute to the electron and 
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FIG. 9: The (mL^/j)3-dependent constraints on from neutron and Mercury EDMs. 



to the quark EDMs, Chromo-EDM s, 3-gluon operator, and 4-fermion operators at the two- 
loop level, and they all induce contributions to the EDMs of the neutron. Thallium and 
Mercury atoms. Our numerical study indicates that the current Thallium EDM bound is 
in general weaker than neutron and Mercury EDM bounds in constraining 0t, and does not 
put any constraint on 0;,. The relative strength of the neutron and Mercury EDM bounds 



TABLE III: Summary of how the CP-violating phases in MSSM are constrained by current EDM 
bounds of neutron, Thalhum, and Mercury atom, for both case I (with hght first-two generations 
of sfermions as in Eq. [3|) and case II (with heavy first- two generations of sfermions as in Eq. S]). 
We refer as "weakly constrained" if the phase can reach 7r/2 with relevant mass scales within about 
1 TeV, and "strongly constrained" otherwise. 





CASE I 


CASE II 


phases 


d-vi 


dn 




d-vi 


dn 






weakly 


weakly 


weakly w. small tan/3 


not 


not 


not 


02 


strongly 


strongly 


strongly 


weakly w. small tan/3 


weakly w. small tan/3 


not 


03 


not 


strongly 


strongly 


not 


weakly w. small tan/3 


weakly 


0e 


weakly 


not 


not 


not 


not 


not 


0u 


not 


weakly 


strongly 


not 


not 


not 


(t>d 


not 


strongly 


strongly 


not 


not 


not 


0M 


not 


not 


not 


not 


not 


not 


4>c 


not 


not 


not 


not 


not 


not 


0s 


not 


not 


not 


not 


not 


not 


0t 


weakly 


weakly 


weakly 


weakly 


weakly 


weakly 


(l>b 


not 


weakly 


weakly 


not 


weakly 


weakly 


4>T 


not 


not 


not 


not 


not 


not 
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on (f)t and 0b depends on tan/3, which drives the relative size of the Yukawa couphngs of the 
bottom and top quarks. As shown in Fig. [HI at small tan/3, the neutron EDM constraint is 
stronger than that of the Mercury EDM , while at large tan/3, the Mercury EDM constraint 
becomes dominant. The (pt and (pf, are rather loosely bounded and can reach 7r/2 for {mL^Fi)3 
of a few hundred GeV. We find that the 0^ phase is not constrained by current EDM bounds 
at all. 

Similarly, we find that the CP violating phases in the second generation of sfermions 0^, 
0c, and 0s are essentially unconstrained by experimental EDM limits. 

C. Correlated constraints on the most strongly constrained phases 02, 03 and 0„^ 

We summarize how each phase is constrained by current EDM bounds in Table IIIII The 
table shows that more than one phase is strongly constrained in case I (light first generation 
sfermions). These phases include 02,3 in the higgsino-gaugino sector, and (j)u,d in the squark 
sector. In the following, we study the EDM constraints on the parameter space defined by 
pairs of such CP violating phases. For simplicity, we tie 0„,d together, but keep 02 and 03 as 
independent. For each of the three combinations of phases (02, 03), (02, 4>u,d), and (03, (l)u,d), 
we perform analysis and determine the region satisfying the combined neutron. Thallium, 
and Mercury EDM bounds at 95% c.L, for tan/3 = 3 and 60, and {mL,R)i = 200, 500 and 
1000 GeV. 

As shown in Fig. (TDl the allowed region for phases grows with the increase in mass scales 
in most cases. The only exception is the 03 with large tan/3 (lower right panel of Fig. [TOl) . 
which is due to the non-monotonic behavior shown in Fig. [H lower panels. 

The correlation between constraints on 03 and (j)u,d (lower panels of Fig. [TO!) is the 
strongest, since both of them are dominantly constrained by the same bound, namely that 
from the Mercury EDM. On the other hand, the correlation between (02, 03) (Fig. [TO] upper 
panels) and (02, (pu,d) (Fig. [ID] middle panels) are rather weak, because 02 is constrained by 
different bounds - neutron and Thallium EDMs. 

As a final illustration of the impact of the EDM constraints, we show in Table IIVI the 
combined bounds (at 95% c.l.) on all the three phases (02, 03, (pu^d) implied by the current 
EDM experimental limits. For purposes of comparison, we also show the impact of (a) 
the reduction in the Mercury EDM bound compared to the previous result {d^^™^^ = 
2 X 10"^'^ecm(95% c. 1.)) ( see Table [VT) : (b) a neutron EDM bound that is 100 times more 
stringent than the present limit, keeping all other EDM limits as at present fIVip : and (c) 

TABLE IV: Summary of the combined bounds at 95% c.l. on three phases (02, 03, 0M,d) for 
tan/3 = 3, 60 and (mi,/j)i = 200, 500, and 1000 GeV, using current experimental limits of neutron. 
Thallium, and Mercury EDMs as in Eq. (1). 



tan/3 


3 


60 


("iL,_R)l 


200 GeV 


500 GeV 


1000 GeV 


200 GeV 


500 GeV 


1000 GeV 


I02| 


< 2.1 X 10-'^ 


< 5.0 X 10-^ 


< 1.5 X 10-2 


< 9.3 X 10-5 


< 2.5 X 10-^ 


< 6.9 X 10-^ 


I03| 


< 2.8 X 10^2 


< 9.7 X 10-3 


< 2.8 X 10-2 


< 3.1 X 10-"^ 


< 4.2 X 10-^ 


< 1.5 X 10-3 


I0«,d| 


< 1.8 X 10-2 


< 6.0 X 10-2 


< 0.17 


< 1.7 X 10-2 


< 5.6 X 10-2 


< 0.21 
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FIG. 10: The correlated constraints on (1)2 and (/)3 (upper panels), (/)„^rf (middle panels) and on (/)3 
and (j)u,d (lower panels), for tan/3 = 3, 60 and {mL^ji)i = 200 GeV (red solid line), 500 GeV (blue 
dashed line), and 1000 GeV (black dotted line)|.gPoints inside the curve satisfy the 95% c.l. of the 
neutron. Thallium, and Mercury EDM bounds. 



a similar improvement in the Thallium limit, keeping the neutron and Mercury limits as at 
present (IVIip . 

The bounds on phases (02, 03, 4>u,d) from the combined 3-phase analysis, as shown in 
Table IIVI are in fact rather close to the bounds from the combined 2-phase analysis, for 
02 shown in the upper and middle panels of Fig. [TUl and for 03 and (pu.d shown in the lower 
panel of Fig. [101 simply because a strong correlation only exists between 03 and (j)u,d, and 
because the inclusion of 02 does not substantially alter this correlation. 

TABLE V: Summary of the combined bounds at 95% c.l. on three phases (02, 03, 4'u.d) for 
tan/3 = 3, 60 and {mL^ji)i = 200, 500, and 1000 GeV, using current experimental limits of neutron. 
Thallium EDMs as in Eq. (1), and previous limits of Mercury EDM djlf'^™"** = 2 x 10~^^ecm (95% 
c. l.)[52]. 



tan/3 


3 


60 




200 GeV 


500 GeV 


1000 GeV 


200 GeV 


500 GeV 


1000 GeV 


I02| 


< 2.1 X 10~3 


< 5.0 X 10-3 


< 1.5 X 10-2 


< 9.2 X 10-5 


< 2.5 X 10-^ 


< 6.9 X 10-4 


I03l 


< 8.7 X 10-3 


< 2.3 X 10-2 


< 6.1 X 10-2 


< 1.9 X 10-3 


< 1.5 X 10-3 


< 3.6 X 10-3 


0M,d 


< 2.3 X 10^2 


< 6.6 X 10-2 


< 0.18 


< 2.8 X 10-2 


< 5.6 X 10-2 


< 0.23 



TABLE VI: Summary of the combined bounds at 95% c.l. on three phases (02, 03, 4>u,d) for 
tan/3 = 3, 60 and {mL,R)i = 200, 500, and 1000 GeV, using current experimental limits of Thallium 
and Mercury EDMs as in Eq. (1), and future possible improved neutron EDM limit d^*^'"'^ = 
2.9 X 10-28ecm (90% c. 1.). 



tan/3 


3 


60 


{mL,R)i 


200 GeV 


500 GeV 


1000 GeV 


200 GeV 


500 GeV 


1000 GeV 


\H 


< 2.1 X 10-3 


< 4.4 X 10-3 


< 1.3 X 10-2 


< 3.6 X 10-5 


< 2.2 X 10-4 


< 5.6 X 10-4 


I03| 


< 7.4 X 10-4 


< 6.1 X 10-3 


< 1.9 X 10-2 


< 2.8 X 10-4 


< 2.8 X 10-4 


< 7.4 X 10-4 


\<Pu,d\ 


< 1.2 X 10-2 


< 3.0 X 10-2 


< 7.8 X 10-2 


< 6.0 X 10-3 


< 2.2 X 10-2 


< 7.6 X 10-2 



TABLE VII: Summary of the combined bounds at 95% c.l. on three phases (02, 03, 0u,(i) for 
tan/3 = 3, 60 and (mL^/j)i = 200, 500, and 1000 GeV, using current experimental limits of neutron 
and Mercury EDMs as in Eq. (1), and future possible improved Thallium EDM limit S^i'^^'^ = 
9.0 X 10-27 ecm (90% c. 1.). 



tan/3 


3 


60 


{mL,R)i 


200 GeV 


500 GeV 


1000 GeV 


200 GeV 


500 GeV 


1000 GeV 


I02| 


< 2.1 X 10-5 


< 5.2 X 10-5 


< 1.7 X 10-4 


< 9.3 X lO-'^ 


< 2.5 X 10-6 


< 7.6 X 10-6 


|03| 


< 2.8 X 10-3 


< 7.8 X 10-3 


< 2.2 X 10-2 


< 3.1 X 10-4 


< 3.8 X 10-4 


< 1.2 X 10-3 


I0u,d| 


< 1.5 X 10-2 


< 4.6 X 10-2 


< 0.13 


< 1.5 X 10-2 


< 5.0 X 10-2 


< 0.19 
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Comparing Table IIVI and |Vl we observe that the most significant impact of the recent 
update on Mercury EDM bound is on 03, while there is no impact on 02 and only small 
impact on (f)u,d- This can be understood by looking at Fig. [1] (middle and lower panels), 
El and [71 The most stringent bounds on 02 are from the Thallium and neutron EDMs. 
The most stringent bound on 03 is from the current Mercury EDM, and the second most 
stringent bound from neutron EDM is many times looser. The most stringent bound on 0^ 
is from the Mercury EDM bound, but the neutron EDM bound is rather close. 

In principle, one might expect that the impact of future neutron EDM experiments would 
be apparent when comparing Table IIVI and IVII Assuming the neutron EDM limit becomes 
100 times tighter than current one, the limits on all phases 02, 03, and (j)u,d, however, would 
only change to be, at most, a few times smaller. This is because as the neutron EDM 
bound becomes the most stringent one for 02, 03, and (f)u,d, strong correlations among all 
of these three phases emerge, making the limits on all of them much larger than when 
the we consider limits on individual phases. This correlation arises from the presence of 
cancellations between various contributions associated with the different phases, a situation 
that has been noted previously in the literature (see, e.g., Ref. [33l|). 

The correlations arising from this "cancellation mechanism" is most easily observed by 
considering the combined constraints on pairs of phases. To illustrate, we show in Fig. 
[TTl the correlated constraints on (02, 03), (02, 4>u,d), and (03, 0«,(i), respectively, using a 
prospective future neutron EDM bound. Comparing with Fig. [TOl we see that with a 
future neutron EDM bound that is 100 times tighter, the neutron EDM would become the 
dominant constraint. In such a case, the cancellation exists among all the phases 02, 03, 
and (f)u,d, leading to a much narrower region for allowed phase values. We also note that 
the orientations of the 95 % C.L. ellipses involving 03 can differ from what appears in Fig. 
[TOl since for the neutron, the constraints arise from the effects of both the quark EDMs and 
chromo-EDMs, in contrast to the situation for Mercury where the chromo-EDMs dominate. 
The relative importance of the dq and dg contributions to dn can change with mL,R, leading 
to changes in the orientation of the ellipses with the value of these mass parameters. 

To asses the impact of the future Thallium EDM experiments, we compare Table HVl with 
IVIII Assuming the future Thallium EDM bound is 100 times tighter than current one, the 
limit on 02 would correspondingly shrink by approximately 100 times, while the limit on 03 
and (j)u,d would not change much, indicating a 100 times tighter Thallium EDM bound would 
not be more stringent than current neutron and Mercury EDM bounds on constraining 03 
and (j)u,d- 



D. Implication of loosely constrained phases on CP-violating and CP-conserving 
phenomenologies . 

As the EDM sensitivities improve, the impact of the most strongly constrained phases 
on other observables will be reduced. On the other hand, the effects of the more loosely 
constrained phases may still be apparent. In the higgsino-gaugino sector, for example, 
the only loosely constrained phase is 0i (for relatively heavy first generation sfermions 
and small to moderate tan/3). This has important consequences in a rather wide range of 
phenomenological contexts, the most significant one being electroweak baryogenesis. A large 
01 can generate a sizeable contribution to the baryon asymmetry via bino-driven electroweak 
baryogenesis scenario, while the other scenario (wino-driven), which depends on 02, is highly 
constrained by current EDM bounds 
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FIG. 11: The correlated constraints on 02 and 03 (upper panels), (/>„^d (middle panels) and on 03 
and 4>u4 (lower panels), for tan/3 = 3, 60 and {'mL,R)i = 200 GeV (red solid line), 500 GeV (blue 
dashed line), and 1000 GeV (black dotted line). Points inside the curve satisfy the 95% c.l. of the 
future neutron EDM bound, and current ThalftSm and Mercury EDM bounds. 



Moreover, this phase is in general present in processes involving neutralinos both on- 
shell and off-shell. Its effect shifts the neutralino mass spectrum and modify the couplings. 
Among CP-conserving quantities, it can lead to order one changes in the production rate, 
decay width, and branching ratios of neutralinos at colliders [8, i9| . It also modifies the relic 
density, as well as direct and indirect detection rate of neutralino dark matter For 
example, as shown in Ref. j22|, the typical variations in the neutralino relic abundance from 
the 01-dependence of couplings is about C(10 — 100%). 

A potentially more direct probe of 0i may be through its impact on CP- violating observ- 
ables at both the LHC and a future linear collider. At the LHC, it contributes to the triple 
product (pi X P2) ■ ps associated with cascade decays of stops |5H7| 

gg iiii, U -> tx% x] ^ x\l^l~ , (5) 

where the pi are the momenta of final state charged particles in the decay chain of the stop. 
This observable is manifestly T-odd. If strong phases are negligible, then it can provide a 
probe of CPV. It is shown in Ref. [6;] that, assuming an order one phase, the signal can 
be detected with 10^ — 10^ identified events. At a future linear collider, it may contribute 
to the triple product {pi+ x pi- ) ■ pe+ associated with neutralino production and subsequent 
leptonic decays (S- I2I 

e+e-^X?X^, Xl^xll^l', (6) 



which is a genuine CP-odd observable. As shown in Ref. [12[, the CP asymmetry in this 
process can reach 10% for some values of the mass and phase parameters. It also contributes 
to another triple-product (pV x Pe+) ■ sV, a T-odd observable that is constructed using the 
transverse polarization sV of the in the neutralino two-body decay 

e^e-^x\xl Xl^r^r^. (7) 



This correlation has been studied in Ref. [13[, where it is shown that the corresponding 
asymmetry can reach values up to 60%. 

Finally, an off-shell neutralino in loop can in principle generate CP-odd observables in B- 
meson decays. However, the neutralino contributions is in most cases subdominant compared 



to other contributions involving gluions and charginos [li 

Among all the loosely constrained phases (^e,/i,T in the slepton sector and 0c,s,t,fe in the 
squark sector, the implication of the third-generation phases 0^ and (f)t,h are most interesting. 
The phase 0^ contributes to the aforementioned T-odd observable (pV x Ve+ ) ■ sV in the process 
outlined in Eq. ([7]) [13]. The phase (j)t,b may generate sizable effect in B- meson physics, even 
though they are flavor- conserving by themselves. The flavor-violation can be either from 
within the SM (CKM quark mixing), or from beyond the SM (off-diagonal elements in 
squark mass matrices). For example, a large 0j can generate sizable deviation from SM for 
S^{ri')Ksi and generate large CP asymmetries Acpip 57) in 6 — )■ 57, through contributions 
involving charginos and charged-Higgs, where the flavor-violation comes from CKM mixing 
matrix [15]. On the other hand, a deviation from SM for S^fiXs can also be generated by 
chirality- flipping LR and RL gluino contributions [3], which can be induced by a flavor- 
violating chirality-conserving mass insertion (5^^ rr)2z and a chirality-flipping CP-violating 
mass insertion associated with the (ph- {^lrrl)™^^^'^'^ = (^ll rr)'23 x uib^Ab — yUtan/3)/m^. 

The phases (j)t,b also generate mixing between CP-even and CP-odd Higgses in MSSM, 
changing their mass spectrum and couplings jl^, which has important consequences for 
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Higgs searches at the Tevatron and LHC [25]. As another example of their impact on CP- 
conserving observables, the large (f)t,b also changes the neutralino annihilation and scattering 
cross section and thus are important for relic density, as well as for direct and indirect 
detection rates for neutralino dark matter l23|. 




IV. CONCLUSIONS 

In the present study, we analyzed the constraints from electric dipole moments on the size 
of CP- violating phases in the MSSM, utilizing the dominant one- and two-loop contributions. 
We introduced the 2LEDM numerical code, interfaced to CPSuperH2.0, that encompasses all 
these contributions. We pointed out that not all CP violating phases in the MSSM are 
constrained to be small by null results from EDM searches. Our results are summarized for 
the ease of the reader in Table IIIII We differentiate there between a case where the first 
generation sfermions are light, and one-loop contributions to EDMs are significant, and one 
where they are heavy and one-loop contributions to EDMs are consequently suppressed. 

We find that in the gaugino sector, the "wino" phase 02 is fairly strongly constrained in 
both cases of light and heavy first generation sfermion masses {mL,R)i, the "gluino" phase 
(^3 is only strongly constrained in the case with light (771l,r)i; and the "bino" phase (pi can 
be arbitrarily large in the case with heavy {mL^R)i, and can be sizable with light (mi^R)i 
and small tan/3. We note that in earlier studies that employed a universality assumption 
(01 = 02 = 03) this difference in EDM sensitivities to the different CPV phases in the 
gauge-Higgs sector of the MSSM was not apparent. 

Turning to sfermion CP violating phases, only the stop CP violating phase 0t (and only 
more weakly 0b) is marginally constrained in the limits of heavy first generation squarks and 
sleptons. For light first generation sfermions, rather stringent constraints arise for 0^ and 
for 0„ from the neutron and Mercury EDM bounds, and weak constraint for 0e arises from 
the Thallium EDM limit. Virtually no constraint exists from EDMs on the CP violating 
phases for second and third generation leptons and for second generation squarks. 

In summary, we showed in this paper that even in the absence of cancellations between 
contributions from different CP violating phases, null results on searches for the permanent 
EDM of the neutron and of atoms put constraints only on selected CP violating phases, 
leaving ample room for a rich phenomenology related to CP violation at colliders, B factories 
and dark matter searches. 
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